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ABSTRACT 

The  synthetic  aperture  radar  principle  has  been  discovered  in  the  early  50th.  Since  then,  a  rapid 
development  took  place  all  over  the  world  and  a  couple  of  air-  and  space-borne  systems  are  operational 
today.  Progress  made  in  technology  and  digital  signal  processing  lead  to  very  flexible  systems  useful  for 
military  and  civilian  applications. 

Radar  has  proved  to  be  valuable  before,  because  of  its  day-and-night  capability  and  the  possibility  to 
penetrate  clouds  and  rain.  Optical  instruments  however  had  great  advantages  in  the  interpretation  of 
depicted  objects.  The  great  wavelength  of  radar  signals  limits  the  achievable  resolution  in  cross  range 
direction  of  real  aperture  radar  systems.  Thus,  imaging  cannot  be  realized  using  static  radar  systems1. 
The  idea  of  SAR  was  to  transmit  pulses  and  store  the  scene  echoes  along  a  synthetic  aperture  (i.e.  the  path 
of  the  SAR  sensor)  and  to  combine  the  echoes  afterwards  by  the  application  of  an  appropriate  focussing 
algorithm.  The  combination  is  carried  out  coherently. 

As  we  will  see,  it  is  quit  easy  to  understand  the  basic  idea  of  SAR.  We  will  show  also  the  hardware  concept 
or  a  SAR  system  and  give  an  idea  for  a  processing  algorithm. 


1.0  GENERAL  FACTS  ABOUT  SAR 

Today,  synthetic  aperture  radar  (SAR)  plays  an  important  role  in  military  ground  surveillance  and  earth 
observation.  Since  the  late  eighties  a  couple  of  SAR-systems  have  been  developed  for  both  space  and 
airborne  operation.  The  underlying  radar  principle  offers  advantages  compared  to  competing  sensors  in 
infrared  or  visible  spectral  area. 

This  chapter  describes  the  field  of  applications  where  SAR  can  be  used  to  gain  valuable  information. 
Physical  conditions  like  propagation  of  electromagnetic  waves  and  scene  reflectivity  affect  the  choice  of 
several  radar  design  parameters. 

1.1  Applications 

SAR  has  got  a  broad  range  of  applications.  For  remote  sensing  a  couple  of  earth  observing  satellites  are 
currently  in  operation,  having  imaging  sensors  working  in  different  spectral  areas.  The  utilisability  of 
optical  sensors  depends  not  only  on  daylight  but  also  on  the  actual  weather  conditions.  Clouds  and  strong 


1  Later  on,  we  will  see  that  high  resolution  radar  imaging  requires  a  relative  motion  between  the  sensor  and  the  object/scene  to 
be  imaged.  It  is  not  relevant,  whether  this  relative  motion  results  from  scene  or  object  motion. 
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rain  are  impenetrable  for  this  wavelength.  Infrared  sensors  which  are  applicable  day  and  night  are  even 
more  sensitive  on  weather  conditions.  Consequently,  radar  sensors  represent  a  completion  of  the  sensor 
collection  for  remote  sensing.  Beyond  the  overall  availability  of  SAR  images  there  are  further  pros  for  the 
utilization  of  radar.  The  coherent  nature  of  SAR  enables  the  user  to  process  images  of  subsequent 
overflights  for  interferometrical  analyses.  Depending  on  the  radar  wavelength  the  radar  signal  will  be 
reflected  by  vegetation  or  the  ground  structure.  With  the  choice  of  a  concrete  centre  frequency  of  the  SAR 
sensor,  the  developer  decides  about  the  appearance  of  the  resulting  radar  images.  Different  combinations 
of  the  transmit  and  receive  polarization  can  also  be  used  for  instance  to  classify  the  kind  of  vegetation. 

Also  in  military  context,  the  availability  of  SAR  is  its  convincing  advantage.  Applications  in  this  area  are 
wide  spread:  global  reconnaissance  is  done  mainly  by  satellite  systems,  aircrafts  and  high  flying 
unmanned  platforms  carry  sensors  for  wide  area  observation  and  miniaturized  SAR  equipment  is  used  for 
integration  into  drones  for  battlefield  surveillance.  The  resolution  of  SAR  systems,  which  indicates  the 
minimal  distance  of  two  small  targets  in  the  scene  which  can  be  separated  in  the  SAR  image,  has  been 
improved  over  the  last  decades  up  to  the  order  of  a  decimetre.  Advanced  classification  algorithms  are  able 
to  identify  military  objects  in  the  scene  which  is  of  great  interest. 

The  exploitation  of  the  relative  motion  between  the  sensor  and  the  scene  is  an  essential  drawback  of  the 
SAR  principle.  While  this  motion  can  be  derived  completely  for  a  fixed  scene  from  the  motion  data  of  the 
SAR  platform,  moving  objects  own  an  additional  motion  component  which  disturbs  the  imaging  process. 
For  future  systems  an  inverse  SAR  mode  (ISAR)  will  be  added  to  image  moving  objects.  No  doubt,  this  is 
an  essential  feature  for  military  application.  In  a  realistic  scenario,  one  sensor  will  perform  the  compound 
task  to  monitor  a  ground  area,  indicate  and  track  moving  objects  and  finally  produce  highly  resolved 
images  of  the  moving  objects  for  target  classification. 

Operation  of  a  bi-  or  multi-static  configuration  offers  some  more  advantages.  The  receiver  system, 
including  expensive  acquisition  electronics,  need  not  to  transmit  any  energy  and  thus  can  be  designed 
hardly  detectable.  Stealth  targets  feature  a  minimized  mono-static  radar  cross  section  (RCS).  Target 
echoes  might  be  considerably  higher  in  a  bi-static  configuration  so  the  probability  of  detection  will  be 
increased. 

To  close  this  discussion  one  related  application  should  be  mentioned.  Ground  based  imaging  of  satellites 
can  be  done  using  the  ISAR  principle.  Optical  sensors  are  inappropriate  for  this  task  because  of  the  great 
distance.  However,  the  satellite  motion  can  be  estimated  and  radar  images  with  high  resolution  can  be 
obtained.  This  is  used  for  damage  assessment  and  space  surveillance. 

1.2  Circumstances 

It  has  been  mentioned  that  the  propagation  and  reflection  of  electromagnetic  waves  strongly  depends  on  its 
frequency.  Figure  1  shows  the  atmospheric  transmission  for  a  wide  frequency  interval.  SAR  systems  for 
wide  area  surveillance  have  to  operate  in  frequency  areas  with  minimal  attenuation.  Beyond  the  Ka-band 
only  short  range  systems  are  realizable.  Advantages  of  systems  with  a  higher  frequency  are  a  reduced 
system  size  and  a  simplified  processing.  X,-  C-  and  L-band  are  very  common  for  airborne  and  spaceborne 
sensors. 

The  wavelength  of  transmitted  signals  plays  an  important  role  also  for  the  reflection  characteristics. 
A  rough  surface  has  completely  different  reflection  behaviour  compared  to  a  smooth  one.  However, 
smoothness  is  a  property  which  depends  on  the  relation  of  the  surface  structure  size  to  the  wavelength. 
Another  feature  of  low  frequency  waves  is  that  foliage  can  be  penetrated  thus  low  frequency  systems  are 
useful  in  some  situations. 
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Figure  1:  Atmospheric  Transmission  of  Electromagnetic  Waves  in  Spectral 
Areas  of  Visible  Light  and  Different  Radar  Bands.  (NASA  1989) 


1.3  Systems  and  Development 

Over  the  years  a  couple  of  systems  have  been  built  and  operated  by  different  nations.  Table  1  shows  some 
systems  to  give  an  impression  of  the  variety. 


Table  1:  Some  Orbital  and  Airborne  SAR  Systems 
Showing  the  Variety  of  Imaging  Radar  Systems 


Shuttle  missions 

SIR-A 

1981 

40  m  x  40  m 

L 

SIR-B 

1984 

25  m  x  17  m 

L 

SIR-C 

1994 

10  m  x30  m 

L,  C,X 

X-SAR/SRTM 

2000 

25  m  x  25  m 

C,X 

Satellite  based 

Lacrosse 

1988 

<  1  m  x  1  m 

X 

ERS-l 

1991 

26  m  x  28  m 

c 

J-ERS-1 

1992 

18mx  18m 

L 

RADARSAT 

1995 

10  m  x  9  m 

C 

ENVISAT 

2002 

25  m  x  25  m 

C 

TerraS  AR-X 

2006 

<  1  m  x  1  m 

X 

Radarsat  II 

2005 

3  m  x  3  m 

C 

SAR-Lupe 

2005 

<  1  m  x  1  m 

X 

IGS-2b 

2008 

30  cm  x  30  cm 

X 

Airborne  SAR 

DOSAR 

1989 

<  1  m  x  1  m 

S,C,X,Ka 

CARABAS-II 

1997 

3  m  x  3  m 

VHF 

PAMIR 

2003 

10  cm  x  10  cm 

X 

Lynx 

1999 

10  cm  x  10  cm 

Ku 

MISAR 

2003 

0.5  m  x  0.5  m 

Ka 

RAMSES 

1994 

10  cm  x  10  cm 

P,L,S,C,X,Ku,Ka,W 

MEMPHIS 

1997 

20  cm  x  20  cm 

Ka,W 

E-SAR 

1994 

1.5  m  x  1.5  m 

P,L,S,C,X 
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2.0  SYNTHETIC  APERTURE  RADAR 
2.1  Principle 

This  section  shows  some  basics  of  the  SAR  principle.  The  presentation  comprises  the  fundamental  facts 
which  will  be  declared  by  descriptive  considerations.  Mathematical  details  are  avoided  for  the  benefit  of 
understanding. 

Assume  a  planar  formation  as  depicted  in  Figure  2  a),  given  in  the  coordinates  x  and  y.  A  radar  sensor  at  a 
known  location  on  the  x-axis  transmits  a  short  pulse2  and  receives  the  echoes  reflected  by  an  object  of  the 
scene.  The  SAR-system  stores  the  received  signals  in  a  two-dimensional  data  array  which  is  parameterised 
in  radar  position  and  echo  signal  delay  which  is  denoted  by  t  in  the  figure.  The  proportional  relation  of  the 
delay  t  and  the  object  distance  r  allows  the  use  of  the  distance  parameter  r  instead  of  the  parameter  t. 
According  to  the  distance  of  all  scene  elements  the  echoes  superpose  each  other  and  result  in  the  recorded 
data  column.  The  data  column  contains  a  range  profile. 


X  x 

a)  b) 


Figure  2:  Two-Dimensional  SAR  Geometry:  a)  data  column  of  a  distributed 
scene  for  a  fixed  antenna,  b)  data  of  point  targets  for  a  moving  antenna. 

Below,  we  consider  a  simplified  scene  as  depicted  in  Figure  2  b).  Three  point  targets3  are  given  at  different 
positions.  The  antenna  moves  in  the  direction  of  the  x-axis.  Since  the  velocity  of  the  radar  signals  is  very 
large  compared  to  the  sensor  velocity,  the  geometry  can  be  assumed  to  be  static  for  one  transmit-receive- 
cycle.  Quasi,  the  antenna  moves  steplike.  A  hyperbolic  range  history  in  the  data  array  results  for  each 


2  Real  systems  use  extended  pulses  which  will  be  compressed  in  the  processing  chain.  For  clarity  we  disregard  this  fact  here. 

3 

A  point  target  is  a  small  object  with  a  reflectivity  which  can  be  assumed  to  be  at  one  discrete  position.  Comer  reflectors  are 
often  used  as  point  targets. 
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reflector  like  it  is  indicated  by  the  curves  in  Figure  2  b).  Points  will  be  generated  in  a  second  data  array  at 
the  positions  of  the  hyperbola  vertex.  If  the  signal  intensity  of  the  individual  echoes  which  results  from  the 
reflectivity  of  the  scene  points  controls  the  brightness  of  the  points  in  the  second  data  array,  an  image  of 
the  scene  results.  For  a  realistic  scene  one  can  imagine,  that  it  is  composed  of  innumerable  point  targets. 
Optimized  algorithms  for  image  generation  focus  the  distributed  energy  of  the  raw  data  for  all  scene 
locations  simultaneously. 

Now,  we  will  discuss  the  achievable  resolution  which  denotes  the  maximal  distance  of  two  identical  point 
targets  which  can  be  separated  in  the  SAR  image.  Since  slant  range  resolution  (resolution  in  y-direction) 
and  the  resolution  in  flight  direction  (x-direction)  depend  on  different  facts  we  discuss  them  one  after  the 
other. 

Using  short  unmodulated  pulses  like  in  this  contribution  leads  to  a  range  resolution  which  is  proportional 
to  the  pulse  duration.  In  the  case  of  rectangular  pulses  the  pulse  duration  is  reciprocal  the  signal 
bandwidth.  A  precise  analysis  shows  that  the  essential  criterion  for  the  achievable  resolution  is  not  the 
duration  but  the  bandwidth  of  the  processed  signal  if  the  possibility  of  pulse  compression  is  exploited. 
Today,  radar  systems  commonly  use  long  pulses  with  large  bandwidth  because  they  offer  high  resolution 
on  one  hand  and  a  convenient  power  ratio. 

The  distance  which  has  been  covered  by  a  reflected  radar  signal  can  be  determined  very  precise  except  for 
a  multiple  of  its  wavelength.  Thus,  the  range  changes  caused  by  changes  of  the  geometry  can  be  measured 
using  periodical  pulse  repetition.  The  accuracy  of  this  measurement  is  better  the  smaller  the  signal 
wavelength  is.  The  thickness  of  the  hyperbolas  in  Figure  2  b)  can  be  chosen  proportional  to  the 
wavelength  to  illustrate  the  influence  of  the  measurement  accuracy.  If  we  analyse  the  hyperbolic  range 
histories  of  two  point  targets  slightly  shifted  in  x-direction,  two  factors  can  be  identified  which  affect  the 
separability  of  the  hyperbolas.  If  the  lines  become  too  thick,  the  hyperbolas  merge  and  cannot  be  divided. 
On  the  other  hand,  limitations  of  the  data  array  in  flight  direction  as  is  depicted  in  Figure  2  b)  by  the 
dashed  box,  cause  a  reduced  separability  of  the  targets.  A  mathematical  analysis  confirms  these  thoughts 
with  the  result  that  the  resolution  in  flight  direction  is  approximately  proportional  to  the  quotient  of  the 
wavelength  and  the  aspect  angle  interval  from  which  the  radar  signals  are  reflected  by  the  scene. 
Two  consequences  should  be  mentioned:  Since  the  observation  angle  cannot  be  arbitrarily  increased, 
the  resolution  in  flight  direction  is  limited4.  Operating  a  SAR  system  in  the  common  stripmap  mode  with  a 
fixed  antenna  look-direction  relating  to  the  flight  direction,  the  antenna  beam  pattern  determines  the 
observation  angle  (equals  the  angle  of  beam).  A  range  independent  resolution  results.  The  observation 
time  of  an  individual  scene  scatterer  (reflecting  element)  increases  proportional  to  the  distance. 

Now  we  extend  the  presentation  to  a  three-dimensional  geometry  to  check  what  happens  in  the  case  of 
imaging  a  special  scene  reflectivity.  Scatterers  of  the  scene  which  will  be  passed  by  the  SAR  sensor  at  the 
same  time  in  the  same  distance  feature  identical  range  histories  and  therefore  superpose  in  one  image 
pixel.  The  interpretation  of  SAR  images  has  to  account  for  this  effect. 

The  design  of  modern  SAR  systems  demands  for  increasing  resolution  capabilities.  The  limitation  which 
is  induced  by  the  operation  in  stripmap  mode  has  become  a  problem.  By  using  a  small  antenna  the  angle 
of  beam  and  so  the  observation  angle  might  be  increased.  Thereby,  the  antenna  gain  would  decrease  and 
the  operational  range  recedes.  Moreover,  a  short  antenna  demands  for  short  distances  between  subsequent 
pulses  which  also  leads  to  range  limitations.  Thus,  the  realisation  of  long  range  systems  is  impossible 
using  this  operational  mode.  The  operation  in  spotlight  mode  solves  this  problem:  While  the  sensor 
platform  passes  the  scene,  the  antenna  beam  will  be  steered  so  that  the  scene  will  be  illuminated  over  a 
large  aspect  angle.  By  this  means  both  demands  can  be  satisfied:  The  operation  of  a  high  gain  antenna  and 


4 


The  theoretical  resolution  which  can  be  reached  observing  the  scene  over  1 80°  aspect  angle  is  a  quarter  of  the  wavelength. 
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the  observation  over  a  large  aspect  angle.  Consequently  most  SAR  systems  currently  under  development 
incorporate  the  spotlight  mode. 

2.2  Imaging  in  Range  Direction 

As  mentioned  before,  modern  radar  systems  transmit  broadband  signals  and  compress  the  received  echoes. 
This  section  describes  imaging  in  one  direction,  concrete  in  range  direction  and  incorporates  pulse 
compression.  Also,  the  presentation  introduces  the  notation  in  spatial  and  wavenumber  domain  which 
turned  out  to  be  practical. 

An  antenna  transmits  the  signal  y,  ( t )  which  will  be  reflected  by  a  point  target  with  reflectivity  a  having  a 
distance  r  from  the  sensor.  The  signal  propagates  with  velocity  c .  The  received  signal  Z(t )  is  composed 
of  the  echo  Jr  ( t )  =  ocst  ( t  -  z)  and  thermal  noise  N(t )  .  Here,  z  —  2r7  c  is  the  signal  propagation  time. 


Z(t)=Sr(t)  +  N(t) 


This  simple  model  neglects  for  propagation  losses.  A  filter  with  impulse  response  h  ( t )  compresses  the 
received  signal.  A  criterion  for  the  efficiency  of  the  combination  of  the  transmit  signal  and  the 
compression  filter  is  the  output  of  the  filter  if  the  transmit  signal  is  compressed.  This  results  in  the  point 
target  response  (also  known  as  point  spread  function) 

p(f)  -  (h  *  st)  ( t ). 


The  asterisk  is  used  as  the  convolution  operator.  Several  aspects  have  an  influence  on  the  choice  of  the 
compression  filter  function.  The  optimisation  of  the  signal-to-noise  ratio  for  t  =  z  is  presented  in  [11]. 

The  result  is  the  matched  filter  which  has  the  impulse  response  hm  ( t )  =  lst  {-t) .  In  this  case,  the  point 
target  response  becomes  the  autocorrelation  function 


Pm 


st{t')st(t  +  t')dt' 


of  the  transmit  signal.  If  5 t(o) ),  p(co)  and  h{co)  denote  the  transmit  signal,  the  point  spread  function 
and  the  transfer  function  of  the  compression  filter,  respectively,  the  filter  function  can  be  expressed  as 


*  **  w 

p  (w)  =h  (u)  st  {u}). 


The  bullet  over  the  signal  symbol  denotes  that  the  function  is  given  in  frequency  domain.  In  the  case  of 
matched  filtering  we  get 


hm 

where 

The  asterisk  denotes  the  conjugate  operator.  Thus,  the  spectrum  of  the  transmit  signal  determines  the  point 
spread  function. 


M  (w)* 


Pm  M  ='?t  St  M  =  st  (u ;) 
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Now,  we  alter  the  presentation  to  spatial  and  wave  number  domain.  Mathematical  expressions  will 
become  more  compact  and  clear.  According  to  the  time  signals  ( t )  and  h  ( t )  we  define 

st(r)  =  st(2rjc)  and  h(r)  =  A(2r/c) 

and  the  point  spread  function  p(r)  =  {h*  st  )(r)  =  cp(2r  /  c)  /  2  . 

Usually,  the  wave  number  kr  of  a  signal,  having  the  wavelength  X  is  determined  by  2k  I X .  However,  in 
the  radar  case  two-way  propagation  has  to  be  taken  into  account,  so  an  equivalent  wavelength  kr  =  4k  /  X 
will  be  used.  In  this  case,  transformation  relations  are  expressed  by 

s  (Ay)  =  r  s{v)e-^rdv 

and 


The  point  spread  function  in  the  case  of  matched  filtering  becomes 

2 


Pm (r)  =  j  |*(  (K 


ejkrrdkr. 


Obviously,  the  spectral  power  density  |i,(£r)|  results  in  a  weighting  in  the  wave  number  domain. 

The  range  imaging  problem  of  a  point  target  in  distance  r'  without  noise  leads  to  the  received  signal 
z(r)  =  ast  ( r  -  r') .  The  output  signal 


zh(r)  ~  (h  *  z)(r)  -  ap(r  -  r') 

of  the  system  for  an  arbitrary  compression  filter  h{r )  .  This  is  a  scaled  and  shifted  version  of  the  point 
spread  function.  Finally,  we  can  extend  the  presentation  to  a  continuous  reflectivity  distribution. 
A  reflectivity  a(t )  returns  echoes  of  the  transmit  signal  st(r )  which  superpose  to  the  signal 

z{r)  =  f  a(/)st(r  —  r/)drf  =  (a  *  $f)(r) 


Application  of  a  compression  filter  with  impulse  response  h(r)  leads  to 

zh(r)  —  (h  *  z)(r)  —  ( h  *  a  *  st}(r)  —  (p  ^  a)(r). 

The  filtered  receive  signal  zh{r )  also  known  as  range  profile  results  as  the  convolution  of  the  scene 

reflectivity  and  the  point  spread  function.  In  the  wave  number  domain  this  expression  is  a  simple 
multiplication  of  the  two  signals. 
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2.3  SAR  Data  Model 

SAR  exploits  the  sensor  motion  to  span  a  large  aperture  to  reach  a  resolution  which  is  much  finder  than 
the  resolution  of  the  real  antenna.  The  data  acquisition  takes  place  in  two  directions.  The  first  direction  is 
range  which  was  discussed  in  the  previous  section.  Repetition  of  the  cycle  Transmit  pulse  -  receive 
echoes’  during  the  sensor  motion  constitutes  the  second  dimension.  Different  parameters  are  authoritative 
for  both  directions:  While  the  velocity  of  light  plays  an  important  role  in  range  direction,  the  sensor 
motion  is  important  in  flight  direction.  Because  of  these  relations,  slow  and  fast-time  are  well  known 
terms  for  the  time  parameters  in  flight  and  range  direction,  respectively.  Following  on,  we  develop  a  signal 
model,  incorporating  the  sensor  motion. 

2.3.1  General  Model 

The  sensor  platform  moves  along  a  path  which  is  described  by  the  vector  b  .  Along  this  path,  an  antenna 
transmits  pulses  and  receives  the  echoes  of  a  scene.  We  suppose,  that  the  platform  retains  its  position 
during  signal  propagation.  Thus,  the  Doppler-effect  can  be  neglected  for  single  pulses.  The  vector 
p  denotes  a  position  of  the  scene,  where  an  echo  will  be  reflected. 


Figure  3:  Geometry  for  General  Model. 


The  signal  model  contains  several  terms: 

•  The  amplitude  of  the  propagated  electromagnetic  wave  decreases  inverse  proportional  to  the 
distance.  After  the  reflection,  this  effect  happens  in  an  equal  way,  thus  it  has  to  be  incorporated  in 
the  received  signal  quadratically. 

•  The  antenna  of  a  SAR  system  has  to  account  for  a  couple  of  different  criterions.  Such  an  antenna 
has  an  antenna  beam  which  leads  to  an  individual  illumination.  This  illumination  of  the  scene  can 
be  described  by  the  function  ga  (p',b)  ,  where  b  is  a  position  of  the  sensor  path.  The  fact,  that  the 

antenna  characteristics  influence  the  signal  twice  has  been  taken  into  account  in  the  choice  of  ga  . 

•  Equivalent  to  the  one-dimensional  presentation,  the  reflectivity  a(p)  denotes  the  ability  of  the 

scene  to  reflect  an  incident  signal.  The  aim  of  the  data  processing  is  to  extract  information  about 
this  reflectivity  from  the  recorded  raw  data.  The  depiction  of  the  reflectivity  is  the  SAR  image  of 
the  scene. 

•  The  transmit  signal  will  be  shifted  according  to  the  distance  of  the  radar  sensor  at  the  current 
position  and  the  target  location. 
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All  partial  echoes  superpose  to  the  scene  echo 

sa&r)  =  JJj  a(p)st{r  -  da(p,b))dp, 

P  ~ 

where  da  (p,b)  =  \p  -  b\  is  the  distance  of  the  two  positions.  The  index  a  marks  the  general  model. 
In  flight  direction,  the  sampling  results  from  the  repetition  of  the  transmit  receive  procedure. 


2.3.2  Model  for  a  Straight  Line 

Most  algorithms  for  SAR  processing  presume  a  straight  flight  path.  Deviations  will  be  handled  as  motion 
errors.  However,  there  exist  also  historical  reasons  for  the  constraint  of  the  straight  flight:  The  idea  of  SAR 
was  to  virtually  create  a  large  aperture  comparable  to  a  linear  array  antenna. 

Given  a  base  location  b0 ,  a  directional  vector  bd  and  a  path  parameter  k  ,  a  straight  plight  path  can  be 
described  by 

b{K)  ~  b0  +  Khj 


The  sensor  approaches  a  point  p  of  the  scene  at  path  position  b(K±  )  if 
(.P-b0j% 


K±  = 


with 


l£rfi 


At  this  position,  the  scene  point  lays  exactly  broadside  the  sensor.  The  distance  between  the  sensor  and  the 
scene  point  can  be  calculated  for  other  path  positions  by 

da(p,b(K))  =  | p-  b(n) |  =  |(p  -  b±)  -  (k  - 


=  da(p,b±)J].  +  (k-«l)2 


l^l2 


(%{p,b±)' 


As  mentioned  before,  the  result  is  a  hyperbolic  range  history.  For  the  special  case  of  a  motion  in 
x-direction,  a  base  location  with  zero  y-component,  a  normalized  direction  vector  as  well  as 

P  =  {px ; Py ;  Pz )  =  {bx  ;0;  bz  Y  and  b0  =  (b0  x  ,0,  b0  z  f ,  some  relations  follow: 

^ X  Px  ^0,x 


d(p,  bx)  =  da(p,b)  =  pp,  1  +  — 

V  TPp 


and 


Pp  =  da{p,bL)  =  yjpy2  +  (pz  -  &0,2)2. 

The  illumination  function  g[p\bx )  results  from  ga  by  evaluation  along  the  straight  path.  The  index  g 
denotes  the  model  for  the  straight  flight. 
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2.3.3  Imaging  of  a  Point  Target 

The  presentation  can  be  reduced  to  point  targets  since  the  received  signal  is  a  superposition  of  their 
echoes.  Insight  can  be  transferred  to  complex  scenes  if  signal  processing  is  also  linear.  The  scene  might 
consist  of  one  point  at  position  PQ=  (/>0  x->Po  y>Po  z)  with  reflectivity  a0  ( p )  =  aS[p^  -  p  ^  ].  The  function 

S{r)  is  the  Dirac  delta.  The  echo  of  the  scene  becomes 

%(£>*,  r)  -  JjJ  ^^a{p)st{r  -  d{p,  bx))dp 

P  " 


A  two-dimensional  compression  can  be  realized  by  applying  a  two-dimensional  matched  filter  which  can 
be  derived  from  the  reference  function.  The  reference  function  can  be  derived  from  the  expected  result  of 
a  point  target  at  the  location  q  having  the  same  coordinates  in  y-  and  z-direction,  but  at  zero  x-position. 
If  g  is  invariant  according  to  translation  in  x-direction,  the  reference  function 


g(q;bx) 

Srcf(bx,r)  =  -  s<(r  -  d{q,  &*)) 

d  [q.  bx) 


(APq^x+Po,*) 

+  Po,*)’ 


St{r  -  di^bx  +  Pa,x)) 


results.  In  the  two-dimensional  wave  number  domain,  the  matched  filtering  can  be  written  as 

Z  m  ( k x  5  kr )  —  $  ( kx  5  kr )  $  ref  ( kx ,  hr )  —  QL  |  S  rej  (  hx ,  hr )  |  C  ^  ^  ’  . 

The  two  closed  bullets  above  the  signal  letters  denote  that  both  parameters  of  the  functions  are  given  in 
wave  number  domain.  Later  we  will  also  work  with  functions  having  one  spatial  and  one  wave  number 
domain  argument.  In  addition,  we  will  then  use  a  bullet  and  a  circle  instead.  The  image  of  the  point  target 
results  after  an  inverse  Fourier  transform: 


Ztn(X)  r)  =  <*Pm(x  PO r) 


with 


Pr 


r)  J~ 2  &  ref  ^ 


The  absolute  value  of  the  two-dimensional  point  spread  function  pm(x,r)  has  a  maximum  at  the  origin 
which  grows  with  the  covered  area  of  the  reference  function  in  wave  number  domain.  The  image  of  the 
considered  point  target  thus  appears  at  position  (/>0x,o)  which  reproduces  the  correct  location  in 

x-direction.  Position  0  in  range  direction  means  here,  that  the  imaged  point  has  the  same  distance  than  the 
reference  point  which  was  assumed  to  produce  the  reference  function.  The  resolution  which  can  be 
achieved  in  flight  direction  depends  on  the  covered  wave  number  domain  area  in  kx .  We  introduce 

Kx  =  max^^-min^),  the  length  of  the  covered  interval.  The  signal  phase  of  the  reference  function 
for  a  mono-frequent  transmit  signal 


S-pef  ^  ) 


cP(pfl,bx  +po,*) 


e  j  ( r  -  d(py  Tbj.  +/J0  )  ) 


can  be  used  to  determine  the  wave  number 
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MM  =  -  Mv  K  +  po,x))  =  -fcr,o-  i — == 

db' 


where  p0  is  the  distance  of  the  analysed  point  to  the  flight  path.  The  antenna  characteristic  can  be 

approximated  by  a  rectangular  function.  If  we  denote  the  angle  between  the  flight  direction  and  the 
direction  to  the  point  target  by  q>  and  its  cosine  by  u ,  this  characteristic  can  be  written  by 

H=lcosW^“*  .The  illumination  function  g[p^^bx  +  p0  x  ]  thus  can  be  converted  to 

IM  <  U9]/p'lp  +  bl 

and 

"y”  =  br{)Ug  ^  kx(bx^  ^  brfiXig  =  2 

limits  the  ‘visible’  area  in  wave  number  domain.  This  rectangular  weighting  in  wave  number  domain  leads 
to  an  si-function  in  spatial  domain,  after  the  inverse  Fourier  transform.  The  Rayleigh  resolution  is 
8X  ~ 71  oug)=  •  Obviously,  the  resolution  is  independent  of  the  point  target  distance. 

The  assumptions  made  here  represent  the  operation  in  stripmap  mode. 


2.3.4  Signal  Model  in  Wave  Number  Domain 

We  have  shown  that  the  covered  area  in  wave  number  domain  plays  an  important  role  for  the  imaging 
capability  of  a  SAR  system.  Moreover,  the  last  section  demonstrated  that  the  compression  can  be  realized 
efficiently  in  the  wave  number  domain  using  simple  multiplications.  This  fact  motivates  to  analyse  the 
data  in  the  two-dimensional  wave  number  domain.  The  Fourier  transformation 


/OG  AOQ 

/  s(bx,r)e-j^T+k^drdbx 

■OG  J  —  GO 


of  the  echo  of  a  point  target  can  be  determined  analytically.  The  Fourier  transformation  in  range  direction 
leads  to 


0* 

s 


(M  hr)  =  r  s(bx,  r)e~jkr 

J  —  OG 


=  Q 


d2(Po>M 


-3krrdr 

st  (, kT)e~jkrd fe’6*) 


The  next  step  converts  the  signal  in  wave  number  domain  also  in  the  direction  of  x.  We  evaluate 

S  (M  fcr)  =  r  8  (bx,kr)e-jk*b*dbx 

J  —  oc 

An  exact  expression  cannot  be  determined.  However,  the  principle  of  stationary  phase  (see  [2]  for 
information  about  the  principle  of  stationary  phase)  enables  us  to  approximate  the  integral.  The  idea  is  to 
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find  one  stationary  point  bx  which  particularly  governs  the  value  of  the  integral.  Other  positions  did  not 
contribute  noticeable  because  of  fast  fluctuations  of  the  phase  term.  The  bx  -dependant  phase  term  of  the 
integrand  is 

*(SiA)  -  -Mfo,,  b,)  -  kxbx 


Figure  4:  Illustration  of  Phase  Terms  to  Interpret  the  Principle  of  Stationary  Phase. 


Figure  4  depicts  the  addends  of  the  phase.  If  kr  is  chosen,  the  contour  lines  of  krd[p^,bx)  are  circles  in 
the  depicted  plane.  The  contour  lines  of  kxbx  are  dashed  lines  in  the  figure.  The  dotted  line  marks  the  path 
of  integration  for  one  concrete  pQ  p .  Along  this  line  we  find  the  stationary  point  where  the  contour  lines 
possess  the  same  distance  in  x-direction.  Stationary  points  fulfil  the  condition 


d<t>(bx) 


dbx 


MpJ  -pm 


sAn)  (Mp0)  -  pm)2 


and  we  find  after  some  steps  one  stationary  point  for 


Mp0)  -  Pcu 


The  line  in  Figure  4  depicts  the  locations  of  stationary  points  for  a  fixed  p0  x  and  variable  p0  . 
The  stationary  phase  for  this  point  results  in 

Mli,))  -  *A(£o) 

k2  - 

=  xr(b*(£n)  -  pm)  -  fcA(Po) 

K'X 

k2  —  k2 

=  —  ,  *  fe(Pn)  -  Po,x)  -  fc*P0,* 

=  —  yjk2  —  k2pQp  —  kxpQjX, 
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Thus,  we  get  the  following  wave  number  domain  model  for  the  echo  signal  of  a  point  target: 


s  ( kx,kr )  =  OTi  st  (kT)- 


d2  ( q 


-kxPo,pe-jkxPo,X 


The  amplitude  term  wx  results  from  the  principle  of  stationary  phase.  However,  this  term  plays  a  minor 
role  in  SAR  processing. 


2.4  SAR  System 

Main  parts  of  a  SAR  system  are  depicted  in  Figure  5.  A  pulse  generation  unit  creates  pulses  with  a 
bandwidth  according  to  the  aspired  range  resolution.  They  will  be  amplified  by  the  sender  and  are 
transferred  to  the  antenna  via  a  circulator.  The  receiver  gets  the  antenna  output  signal  (echoes  of  the  scene) 
amplifies  them  to  an  appropriate  level  and  applies  a  bandpass  filter.  After  the  demodulation  and  A/D 
conversion  of  the  signals  the  SAR  processor  starts  to  calculate  the  SAR  image.  Additional  motion 
information  will  be  provided  by  a  motion  measurement  system.  A  radar  control  unit  arranges  the  operation 
sequence,  particularly  the  time  schedule. 


Figure  5:  General  Structure  of  a  SAR  System. 


2.5  Processing  of  Radar  Data 

Over  the  years,  a  couple  of  processing  algorithms  have  been  developed.  At  first,  separated  processing 
approaches  have  been  implemented,  which  sequentially  solved  the  problems  of  range  and  azimuth 
(in  flight  direction)  compression.  These  methods  are  well  documented  in  the  literature  as  range  Doppler 
processing.  With  the  increasing  performance  of  newly  developed  systems,  these  processing  algorithms 
reached  their  limitations  which  originate  from  implied  approximations.  Latest  developments  are  useful  to 
process  data  up  to  any  resolution.  However,  the  efficiency  has  to  be  bought  by  extensive  hardware 
investment. 

Following  on,  we  will  present  one  processing  strategy.  The  signal  model  of  a  point  target  in  wave 
number  domain  has  been  derived  in  section  2.3.4.  Analysing  the  phase  of  that  signal,  we  find  that  it  is 
non-linear  in  the  wave  number  coordinates.  The  idea  of  the  range  migration  algorithm  (also  known  as 
co-k  -algorithm)  is  to  introduce  a  new  coordinate  in  the  wave  number  domain  so  that  the  signal  gets  a 
linear  phase  in  the  new  coordinate  system.  The  substitution 

kP  =  A2  -k2 
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solves  this  problem.  An  inverse  Fourier  transform  leads  to  the  SAR  image  of  the  point  target. 
(The  amplitude  term  is  the  wave  number  representation  of  the  point  spread  function  and  the  linear  phase 
term  provides  the  correct  location  after  the  transformation.)  Within  the  implementation  the  substitution 
enables  a  resampling  in  the  new  coordinate  k p  .  The  new  samples  will  be  calculated  using  the  so-called 
Stolt-interpolation. 

Further  information  about  SAR  data  processing  can  be  found  in  a  couple  of  books  which  are  given  in  the 
literature  list. 
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